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Abstract: 

An entirely quantum mechanical approach to diagonalize hermitean matrices has been pre- 
sented recently. Here, the genuinely quantum mechanical approach is considered in detail 
for (2 x 2) matrices. The method is based on the measurement of quantum mechanical ob- 
servables which provides the computational resource. In brief, quantum mechanics is able to 
directly address and output eigenvalues of hermitean matrices. The simple low-dimensional 
case allows one to illustrate the conceptual features of the general method which applies to 
(N x N) hermitean matrices. (Fortschr. Phys. 51, 248 (2003)) 

PACS: 03.67.-a, 03.65Sq 
Introduction 

A new attitude towards quantum theory has emerged in recent years. The focus is no longer 
on attempts to come to terms with counter-intuitive quantum features but to capitalize on 
them. In this way, surprising methods have been uncovered to solve specific problems by 
means which have no classical equivalent: quantum cryptography, for example, allows one 
to establish secure keys for secret transmission of information [1] ; entanglement [2] is used 
as a tool to set up powerful quantum algorithms which do factor large integers much more 
efficiently than any presently known classical algorithm [3]. Throughout, these techniques 
make use of the measurement of quantum mechanical observables as an unquestioned tool. 
This is also true for many (but not all [4]) proposals of quantum error correction schemes 
[3, 5], required to let a potential algorithm run. 

The purpose of the present contribution is to study the simplest situation in which a 
quantum mechanical measurement, i.e. the bare 'projection' [2], "does" the computation. 
The computational task is to determine the eigenvalues of hermitean (2 x 2) matrices by 
quantum means alone. Although the answer to this problem can be given analytically, 
it is useful to discuss this particular case since there is no conceptual difference between 
diagonalizing (2 x 2) or (N x N) hermitean matrices along these lines [6]. 

Before turning to the explicit example, consider briefly the traditional view on quantum 
mechanical measurements: a measurement is thought to confirm or reveal some information 
about the state of the system. The measured observable A is assumed to be known entirely, 
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including in-principle-knowledge of its eigenstates and eigenvalues. Further, the observable 
also defines the scope of the possible results of a measurement since both the only allowed 
outcomes are its eigenvalues and, directly after the measurement, the system necessarily 
resides in the corresponding eigenstate. 

In the context of quantum diagonalization, however, the crucial idea is to learn something 
about the measured observable — not about the state of the system. Why is there scope for 
information gain at all? What can one learn from a quantum mechanical experiment if both 
the measured observable and the state are known? 

It is essential to realize that the input required to measure an observable A and the 
output of an experiment, in which A is actually measured, are not identical. In fact, it 
is possible to construct an apparatus which measures an observable A without explicitly 
knowing its eigenvalues. Then, a measurement provides partial — but explicit — information 
about the spectral properties of the observable A: it delivers one of its eigenvalues. As the 
eigenvalues of A are not known explicitly before the measurement, information is indeed 
gained by measuring A. This is the idea which underlies quantum diagonalization. 

The quantum diagonalization of a hermitean matrix is achieved in five steps: (1) express 
the matrix in a standard form] (2) associate a quantum mechanical observable with it; (3) 
identify an apparatus capable of measuring the observable; (4) measure the observable — this 
provides the eigenvalues of the matrix; (5) determine its eigenstates. The next section gives 
the details for (2 x 2) matrices. Then, the generalization to (N x N) matrices is briefly 
summarized. 

Quantum diagonalization of hermitean (2 x 2) matrices 

Suppose you want to determine the eigenvalues A± (and, subsequently, the eigenvectors 
\A±)) of the general hermitean (2 x 2) matrix 

a,j€lR, 0e€. (1) 
Obviously, this problem is easily solved analytically. The eigenvalues read 

A± = \ (a + 7 ± V(a-7) 2 + W*) , (2) 

and one can also give explicit expressions for the eigenvectors of the matrix A. The five- 
step procedure of quantum diagonalization is now applied to A; it will, of course, reproduce 
the result (2). However, all conceptual points, which also apply to the technically more 
cumbersome case of (N x N) matrices are conveniently illustrated by this simple example. 

1. Standard form of A: Any hermitean (2 x 2) matrix A can be written as a unique linear 
combination of the three Pauli matrices Oj = Oj,j = 1,2,3, , and the unit matrix 
co = h, 

A = (a cr + a • a) , 

where 

a = ^(a + 7), a x = 1(0 + /?), 



a 0* 
7 



a j = -Tr [Aa 3 ]eM, 



(3) 



a 2 = ^(/3*-/3), a 3 = i(a-7). (4) 
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From a general point of view, this corresponds to an expansion of A in multipole op- 
erators (cf. below). 



2. Identification of an observable: The most general Hamiltonian of a spin- 1/2 in a homo- 
geneous magnetic field B is linear in the components of the spin S = ha/2. Therefore, 
any matrix A has an interpretation as a Hamiltonian operator of a quantum spin sub- 
jected to an appropriately chosen magnetic field, 

_2 

A = a\ 2 - .g^sBo • S = H A (S) , a = a , B = -a . (5) 

The part ah shifts the energy globally by a fixed amount. Let E± denote the eigenval- 
ues of the second part of this expression, — g^B • S = H^(S); then, the eigenvalues 
A± of the matrix A are given by 

A± = a + E± . (6) 

Consequently, the door is now open to determine the eigenvalues of A experimentally, 
i.e. through measuring the eigenvalues E± of the Hamiltonian H^(S). In the next step 
it is shown how to devise an apparatus which measures this operator. 



3. Setting up a measuring device: The apparatus is required to measure the eigenvalues 
of the operator H^(S). In the case of a (2 x 2) matrix this is just a familiar Stern- 
Gerlach apparatus, appropriately oriented in space. However, as the method will be 
applied to (N x N) hermitean matrices later on, it is important to go through the 
construction of the measuring device in detail. 

Consider the spatially varying Hamiltonian 

H°(r, S) = ah - OMsB(r) • S = al 2 + H° (r, S) , (7) 

which describes the interaction of a spin with an inhomogeneous magnetic field B(r). 
In order that H°(r, S) measure the observable H^(S), the field needs to satisfy the 
conditions 

B(0) = Bo and V • B(r) = V x B(r) = . (8) 

Then, at the centre of the apparatus, the operator H^(r, S) coincides with the ob- 
servable to be measured, H^(0,S) = H^(S), and the magnetic field is consistent with 
Maxwell's equations. Consider the field [7] 

B(r) = (l + k-r)B + (B .r)k, (9) 

which is consistent with Eqs. (8) if the vector k is perpendicular to B . Diagonalizc 
the operator H^(r, S) to first order of V|B|/|B|. This leads to r-dependent eigenvalues 

E ± (r) = ±J(l + kT) J B , (10) 
which imply the existence of a state-dependent force 

F±(r) = -V£7±(r) = T^S k. (11) 
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Consequently, particles on their way through the apparatus will be deflected determin- 
istically from a straight line once the projection to an eigcnstate \E±) has occurred. 
In this way, the eigenvalues of the observable can be accessed experimentally. 

4. Determination of the eigenvalues: If a measurement of the operator H^(S) is performed 
on the state p = h/2, then it is thrown with probability 1/2 into one of the eigenstates 
\E±) with density matrix p±, 

app(H°) : p = l 2 /2 ^ (E±; p± ) , p± = Tr [pp±] = ± . (12) 

Repeating the measurement a few times, both eigenvalues will have been found soon. 
The probability not to obtain one of the two values after No identical runs of the 
experiment equals 1/2 N °, and hence goes to zero exponentially with the number of 
runs. Subsequently, the sought-after eigenvalues A± of the matrix A are known due 
to the relation (6), and the major step in the diagonalization of the matrix has been 
achieved in a quantum way. 

5. Determination of the eigenstates: Once both eigenvalues A± are known, it is straight- 
forward to to determine the associated eigenstates by a simple calculation. Optionally, 
one continues along an experimental line. One exploits the fact that the apparatus 
app(H^) prepares an ensemble of eigenstates of A with density matrix p + if the other 
subbeam (containing p_) is blocked — and vice versa. The Bloch representation of the 
density matrix p + , say, can be parameterized by expectation values, 

P+ = \ (°o + • cr) , (<7j-) + = Tr [p + o-j] = (E + \aAE + ) . (13) 

Hence, the components of the vector (er)+ (and therefore p + ) are easily obtained by 
means of a second, appropriately oriented Stcrn-Gerlach apparatus, which amounts to 
a reconstruction of the density matrix p + . 

Quantum diagonalization of hermitean (N x N) matrices 

Five steps are necessary to diagonalize a hermitean (N x N) matrix A by quantum means. 

1 . Standard form of A: Write the hermitean (NxN) matrix A as a combination of linearly 
independent hermitean multipole operators T v , v = 0, . . . , TV 2 — 1, 

A = a » J v . a » = Tr Tr I AT ^] EM - ( 14 ) 

Multipole operators Tj U2 ...j a act in a Hilbert space H s of dimension N = 2s + 1 
which carries an irreducible representation of the group SU{2) with the spin com- 
ponents (Si,S2,S3) as generators. They are defined as the symmetrized products 
SjjSjv, • • -Sj a ,ji = 1,2,3, and a = 0, 1, ... ,2s, after subtracting off the trace, except 
for T = T(°) = I, the (N x N) unit matrix. The index a labels (2s + 1) classes with 
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(2a + 1) elements each, transforming among themselves under rotations. Explicitly, 
the lowest multipoles read 

tW — Q . t(2) _lfC.C. I c c \ jjUj c C 

1 j — J J i ji32 ~ 2 ^ n 32 ■ 3 32 J 3l) g J Jl J j2 ■ \ LO ) 

For the sake of brevity, a collective index v = (a;ji,... ,jk) has been used in (14), 
taking on the values v = 0,1, . . . , N 2 — 1. The N 2 self-adjoint multipole operators 
T v = Tj, form a basis in the space of hermitean operators acting on the TV-dimensional 
Hilbert space H s [7]. Two multipoles are orthogonal with respect to a scalar product 
defined as the trace of their product: (1/N) Tr [T^T^] = 8 vv i . 



2. Identification of an observable: Interpret the matrix A as an observable for a single 
quantum spin S with quantum number s = (N — l)/2, 

N 2 -l 

A = a -"U S ) = M S ) . (16) 

using the expression of the multipoles T„(S) in terms of the components of a spin. 
Since the multipoles are expressed explicitly as a function of the spin components not 
exceeding the power 2s, it is justified to consider them and, a fortiori, the quantity 
H J 4(S) as an observable for a spin s. 



3. Setting up a measuring device: Swift and Wright [7] have shown how to devise, in prin- 
ciple, a generalized Stern-Gerlach apparatus which measures any observable H^S) — 
just as a traditional Stern-Gerlach apparatus measures the spin component n • S along 
the direction n. The construction requires that arbitrary static electric and magnetic 
fields, consistent with Maxwell's equations, can be created in the laboratory. To con- 
struct an apparatus app(H y i) means to identify a spin Hamiltonian H(r, S) which splits 
an incoming beam of particles with spin s into subbeams corresponding to the eigen- 
values A n . The most general Hamiltonian acting on the Hilbert space Ti. s of a spin s 
reads 

N 2 -l 

H(r,S) = MOT,, (17) 

with traceless multipoles (except for v = 0), and totally symmetric expansion coef- 
ficients $„(r) (= &j\] 2 ...j k ( r ))- These functions, which vary in space, generalize the 
magnetic field B(r) in (17). Tunc the corresponding electric and magnetic fields in 
such a way that the coefficients <!>„ (r) and their first derivatives with respect to some 
spatial direction, r 1; say, satisfy 



3v( r )Uo = a ^ and 



dr\ 



(18) 



r=0 



As shown explicitly in [7] , this is always possible with realistic fields satisfying Maxwell's 
equations. Then, the Hamiltonian in (17) has two important properties: 

(i) At the origin, r = 0, it coincides with the matrix H(0, S) = H^S). 
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(ii) Suppose that a beam of particles with spin s enters the generalized Stern-Gerlach 
apparatus app(H J 4) just described. At its center, particles in an eigenstate \A n ), 
say, will experience a force in the r\ direction given (up to second order in distance 
from the center) by 

d(A n \H{v,S)\A n ) 
^i(r)Uo = Wi 

Consequently, particles with a spin residing in different cigcnstates \A n ) of the opera- 
tor Ha will be separated spatially by this apparatus. From a conceptual point of view, 
the procedure is equivalent to the method outlined above for a spin 1/2. 



= -A n , n= l,...,2s + l. (19) 



4. Determination of the eigenvalues: Once the apparatus app(H J 4) has been set up, one 
needs to carry out measurements on (particles carrying) a spin s prepared in the 
homogeneous mixture p = \n/N. The output of each individual measurement will be 
one of the eigenvalues A n of the matrix A, according to the 'projection postulate:' 

v 1 
app(A) : p=\ N /N (A n ; Pn ), p n = Tr [pp n ] = — . (20) 

In words, the action of the apparatus is, with probability p n = 1/N, to throw the 
system prepared in state p into an eigenstate p n = \A n )(A n \ of the observable A; the 
outcome of the measurement is given by the associated eigenvalue A n . 

After sufficiently many repetitions, all eigenvalues will be known, although the outcome 
of an individual measurement cannot be predicted due to the probabilistic character 
of quantum mechanics. It is necessary to repeat the experiment a number of times 
until all values A n have been obtained. Since the spin s has been prepared initially in 
a homogeneous mixture, p = \n/N, the (2s + 1) possible outcomes occur with equal 
probability. The probability not to have obtained one specific value A n after No ^> N 
measurements equals (2s/(2s + 1))^° ~ exp[— Nq/2s], decreasing exponentially with 
N . 

5. Determination of the cigcnstates: As before it is possible to either calculate the eigen- 
states \A n ) of the matrix A on the basis of the known eigenvalues, or to determine 
them experimentally by methods of state reconstruction (see [8] for details). 

Summary and Outlook 

As a result of the five steps just described, a hermitean (N x N) matrix A has been diagonal- 
ized without calculating the zeroes of its characteristic polynomial by traditional means. The 
fourth step solves the hard part of the eigenvalue problem since it provides the eigenvalues 
A n of the matrix A — information which cannot be obtained in closed form if TV > 5. One 
might best describe the measuring device app(Hyi) as a 'special purpose machine' which is 
based on the 'collapse of the wave function' as computational resource. However, one could 
avoid to use the notion of 'collapse' or 'projection' by characterizing the process indirectly 
using the concept of 'repeatable measurements' described in [9]. 

By construction, the quantum mechanical diagonalization is not based on the representa- 
tion of a mathematical equation in terms of a physical system which then would 'simulate' it. 
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Similarly, no 'software program' runs which would implement an diagonalization algorithm. 
Therefore, the method resembles neither an analog nor a digital calculation. 

t is worthwhile to point out that the quantum mechanical approach to the diagonalization 
of hermitean matrices is based on the assumption that the behaviour of a spin s is described 
correctly by non-relativistic quantum mechanics. Note, finally, that quantum diagonalization 
does not depend on a particular interpretation of quantum mechanics. 
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